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Abstract 

We investigate the relationship between the covariant and the three-dimensional 
(equal-time) formulations of quantum kinetic theory. We show that the three- 
dimensional approach can be obtained as the energy average of the covariant for- 
. mulation. We illustrate this statement in scalar and spinor QED. For scalar QED 

we derive Lorentz covariant transport and constraint equations directly from the 
Klein-Gordon equation rather than through the previously used Feshbach-Villars 
representation. We then discuss the relation to their equal-time version. We make 
a semiclassical expansion in h and obtain Vlasov-type transport equations for the 
zeroth and first order of the scalar Wigner operator. We then consider pair pro- 



X 



duction in a spatially homogeneous but time-dependent electric field and show that 
the pair density is derived much more easily via the energy averaging method than 
in the equal-time representation. Proceeding to spinor QED, we derive the covari- 
ant version of the equal-time equation derived by Bialynicki-Birula et al. We show 
that it must be supplemented by another self-adjoint equation to obtain a complete 
description of the covariant spinor Wigner operator. After spinor decomposition 
and energy average we study the classical limit of the resulting three-dimensional 
kinetic equations. There are only two independent spinor components in this limit, 
the mass density and the spin density, and we derive also their covariant equations 
of motion. We then show that the equal-time kinetic equation provides a complete 
description only for constant external electromagnetic fields, but is in general incom- 
plete. It must be supplemented by additional constraints which we derive explicitly 
from the covariant formulation. 
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1 Introduction 



Theoretic evidence from lattice simulations Jl| of QCD for a phase transition from a 
hadron gas to a quark-gluon plasma at high temperature has prompted experimental 
efforts to create this new phase of matter in the laboratory during the early stages of 
ultrarelativistic heavy ion collisions p. However, because of the estimated very short 
life time of the collision zone, the highly excited quark-gluon system may spend a con- 
siderable fraction of its life in a non-thermalized, pre-equilibrium state. The dynamical 
tool to treat dissipative processes in heavy ion collisions and the approach to local ther- 
mal equilibrium is in principle non-equilibrium quantum transport theory. Since QCD, 
the theory describing the interactions between quarks and gluons, is a gauge theory, the 
kinetic equations should be gauge covariant. A relativistic and gauge covariant kinetic 
theory for quarks and gluons has been derived both in a classical framework j|, |5| and 
as a quantum kinetic theory || [7| based on Wigner operators defined in 8-dimensional 
phase space ||. Some preliminary applications to the quark-gluon plasma, such as lin- 
ear color response, color correlations and collective plasma oscillations || [X0|] have been 
discussed in this framework using a semiclassical expansion of the quantum transport 
theory. Apart from its relativistic and gauge covariance, an important aspect of the four- 
dimensional approach to transport theory is that the complex kinetic equation can be split 
up into a constraint and a transport equation, where the former is a quantum extension 
of the classical mass-shell condition, and the latter is a covariant generalization of the 
Vlasov-Boltzmann equation. The complementarity of these two ingredients is essential 
for a physical understanding of quantum kinetic theory |LT] . 

As a first step towards a full kinetic treatment of the quark-gluon plasma, a detailed 
study of relativistic transport theory for electrodynamics was performed over the years 
(see H for references). The absence of the complications arising from the non-abelian 
character of QCD allows for a deeper insight into the structure of the kinetic theory 
itself, and it is also easier to find useful applications in this simpler case. Vasak, Gyu- 
lassy and Elze |13| (to be quoted as VGE) discussed the relativistic quantum transport 



theory for spinor electrodynamics and gave a very lucid account of the covariant spinor 
decomposition for the Wigner operator of Dirac particles (see also ]n| ) . 

In classical transport theory, all the physical currents are connected with the distri- 
bution function /. The quantum mechanical analogue of / is the Wigner function, which 
is the ensemble expectation of the Wigner operator. In the process of performing the 
ensemble average of the kinetic equation for the Wigner operator, one encounters the 
two-body correlation function (F^W), where is the gauge field strength. In the 
general case, the two-body correlation function again depends on higher order correlation 
functions; this generates the so called BBGKY hierarchy ||. A popular method to obtain 
a closed kinetic equation for the Wigner function W = (W) is to use the Hartree approxi- 
mation where the gauge field is considered as a mean field F^, leading to the replacement 
(F^W) = F^CW). In this approximation the BBGKY chain is truncated at the one 
body level, and the kinetic equation for the Wigner function has the same form as that 
for the Wigner operator. This version of the mean field approximation, which treats the 
particle fields quantum mechanically, but uses the classical approximation for the gauge 
fields, is widely used in electrodynamic transport theory [13|, [14], . It is appropriate for 
strong but slowly varying electromagnetic fields. An interesting example is the Schwinger 
process [fill fL7f for pair production in external electromagnetic fields. Since this approx- 
imation is sufficient for our purpose, we will also restrict ourselves to the investigation of 
transport theory for particles interacting with an external field. 

In a recent paper Bialynicki-Birula, Gornicki and Rafelski [fyj (to be quoted as BGR) 
proposed an equal-time transport equation for spinor electrodynamics. Unlike the covari- 
ant theory where the Wigner operator is defined as the four-dimensional Fourier transform 
of the gauge covariant density matrix <§>±{x,y) (see below), 



they introduced an equal-time correlation function <3> 3 (x, y) = $(x, y,t) and considered 
its spatial Fourier transform which depends only on three momentum coordinates: 




(1) 




(2) 
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From this definition and the Dirac equation they obtained directly a kinetic equation 
for W 3 . They showed that the Wigner function (W 3 (x,p)) is a direct analogue of the 
classical distribution function /(x, p, t) and that it provides a systematic way of studying 
the phase-space dynamics of QED in the semiclassical limit; each spinor component of 
(W3) corresponds to a definite physical distribution function. Of course, Wz(x, p) is not 
manifestly Lorentz covariant. 

Recently this equal-time method has been extended ]TJ| to scalar electrodynamics in 
the Feshbach-Villars representation. The ensuing applications of this so-called "three- 
dimensional" transport theory to the problem of nonperturbative pair creation in both 
spinor |14| and scalar (|Tj| [0| QED indicated that some quantum problems may be treated 



more easily in phase space in terms of the Wigner operator than by conventional field 
theoretic methods. We will show, however, that the published derivation of the equal-time 
transport theory is incomplete. A complete treatment leads to additional constraints on 
the Wigner function which for the specific situations treated by previous authors actually 
simplifies its structure; in general, however, it leads to complications whose effects cannot 
be neglected. 

We show in this paper that the three-dimensional transport theory can also be derived 
by a different method which does not suffer from this incompleteness. In particular, we 
discuss the relationship between the covariant and three-dimensional approaches. From 
the comparison of Eqs. ([I]) and (|2]) it is easily seen that there exists a general connection 
between the four- and three-dimensional Wigner operators: 

W 3 (x,p) = J dp Wt(x,p). (3) 

This relation shows that the three-dimensional transport theory can be obtained by an 
energy average of its covariant version. The main formal difference between the energy 
averaging method used in this paper and the direct equal-time derivation of BGR is that 
our approach starts with an explicitly relativistic formulation. This covariant formula- 
tion is shown to be complete, i.e. the equation of motion for the Wigner operator is 



equivalent to the Dirac equation for the field operators. This completeness is preserved 
during the reduction from the four- to the three-dimensional version via energy averag- 
ing. The covariant formulation followed by energy averaging also allows us to obtain a 
three-dimensional transport equation for the scalar Wigner operator directly from the 
Klein-Gordon equation, rather than via the Feshbach-Villars 2x2 matrix formulation. 

We proceed as follows. In section 2, we discuss the transport theory for spinless charged 
particles. We first set up the covariant constraint and transport equations, which are the 
basis of the reduction to the three-dimensional version. Next we consider the semiclassical 
expansion in H and the classical limit. Then we average the covariant equations over 
the energy to obtain the three-dimensional equations. With these we investigate the pair 
creation of charged scalar particles in a spatially homogeneous but time-dependent electric 



field. This problem was studied before in the Feshbach-Villars matrix representation |19 



but our procedure turns out to be much simpler because of the scalar nature of our 
Wigner function. In Section 3, we study the kinetic theory for Dirac fermions. We derive 
a complete set of two selfadjoint covariant equations, one of which corresponds to the 
four-dimensional version of the BGR equation. We give a full mapping of the spinor 
components from the four- to the three-dimensional representation. The two selfadjoint 
covariant transport equations can be combined into a single complex equation, the VGE 
equation, which is then subjected to the energy averaging procedure. We discuss the 
classical limit of the resulting three-dimensional equations and compare it with the results 



20fl from the equal-time method. We then show how to get a complete set of three- 



dimensional kinetic equations in the general, fully quantum mechanical case; these contain 
as a subset the various spinor components of the BGR equation, but also a set of additional 
constraints. We discuss these results and summarize them in our conclusions. 
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2 Scalar electrodynamics 

2.1 Relativist ic covariant kinetic equations 

The abelian gauge theory of scalar particles with mass m and charge e is defined through 
the Lagrangian density 

C = (0 M - ieAjtfid* + ieA")<j> - m 2 f - \F^F^ . (4) 
It leads to the Klein-Gordon equation for the scalar field operator (ft, 

"(fy + ieA M (ar))(9^ + ieA'*(ar)) + m 2 ] 0(ar) = , (5) 
the adjoint equation for 0', and to Maxwell's equations for the field strength tensor 

«V^(x)=/(x), 9 M F^(x) = 0, (6) 

where F^ u {x) = \e [1V(XJ F pu is the dual field tensor and j u (x) is the charge current of the 
scalar particles. 

The relativistically covariant Wigner operator is the Fourier transform Eq. (|1|) of the 
gauge covariant scalar field correlation function 

$ 4 ( X) y ) = <f)(x)e- y - vt/2 ■ e y - v/2 <p\x) . (7) 

The covariant derivatives 'D^x) = <9 M + ieA^x) and T>Ux) = <9 M — ieA^x) in this ex- 
pression ensure the gauge covariance (actually: gauge independence in the case of a U(l) 
gauge theory) of the density matrix. Following |J and replacing the covariant derivatives 
here by a line integral of the gauge field A^(x) along a straight line between the space-time 
points x — | and x + |, we recover the more familiar form [fL6[ 



i 



$ 4 (x, y) = (f> (x + |y) exp ie / * c?s A(x + sy)-y 0^ (x — \y} . (8) 

L "2 J 

The resulting scalar Wigner operator is selfadjoint: 

Wl(x,p) = W 4 (x,p) . (9) 
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Substituting Eq. ([?]) into Eq. ([!]) and integrating over y we see that the Wigner operator 
is just the density of particles at space-time point x with kinetic momentum p, 



W±(x,p) = 4>{x) 5 i (p - n(x)) 4>\x) 



(10) 



where fr(x) = — DUx)) is the kinetic momentum operator. Eq. (|1(]) motivates 

a statistical interpretation of the Wigner function as a generalized phase-space density 
and generates simple relations between the Wigner operator and all physical space-time 
density operators. For example, the operators for the charge current density and the 
energy momentum tensor are given by 



3A X ) 



ie(ct) ] (x)d ll (i)(x) - (d^ (x))(f)(x)J 
ejd^W^p), 

\([d^{x)){d v ct>{x)) - (d u ^(x))(d,<P(x 
d A pp^p u W 4 (x,p) . 



(11) 

The equations of motion for the Wigner operator are a direct consequence of the field 
equations ((^). We calculate the second-order derivatives of the correlator $4(2;, y) with 
respect to x and y, 

1 

-lie [[ ds (I + s) y v F^(x + sy) (±c£ + S») $ 4 (z, y) 
1 

\ ds (I + s ) y- F ( x + s v) 



-Ue 



$ 4 (x,y) 



-He / 2 (1 + s) y-j(x + sy) $ 4 (a;, y) 



(12) 

where we have employed the Klein-Gordon and Maxwell equations in the first and last 
terms on the right-hand side. Eq. flT2|) can be rewritten in a compact way, 



\d»d^ - rpn M + m 2 - iWD J $4(2, y) = o 



(13) 
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where we denned the two Lorentz covariant operators D and II by 

D»(x, y) = d% + xe f\ ds y v F v »{x + ay) , 

"2 

Wfay) = i(d% + ie J\ds sy u F utl (x + sy^j . (14) 

Performing the Fourier transform with respect to y, we obtain the following exact 
quadratic kinetic equation for the scalar Wigner operator: 

'itfDtDp - ITTL, + m 2 - ihWD^ W 4 (x, p) = . (15) 

The operators D and II in this equation are now defined in phase space and can be 
obtained from Eqs. ( |I~4T) by the replacements y^ — > — id£ and — > —ip^: 

Dp(x, p) = dp-e J x ds F^ u (x - ihsd p ) d v p , 

n M (a;,p) = Pn — ieH ^ ds s F^(x — ihsd p ) d% . (16) 

"2 

In Eqs. QT5D and ( JT6| ) we reinstated the ^-dependence explicitly in order to be able to 
discuss the semiclassical expansion in the next subsection. (The speed of light c is still 
omitted.) Obviously, the operators and 11^ are gauge covariant extensions of the 
partial derivative <9 M and the momentum p^, respectively. They both are self-adjoint: 
namely, D^(x,p) = D^(x,p) and 11+ (a;, p) = Yl^(x,p). 

Since the scalar Wigner operator is self-adjoint, too, the real and imaginary parts of 
the complex equation (|15| ) have to vanish separately: 

(\h 2 D»D^-WIl fl + m 2 )w 4 (x,p) = } (17) 
hW t D li W A (x,p) = 0. (18) 

The first equation is usually called constraint equation since it is obviously a generaliza- 
tion of the classical mass-shell condition p 2 — m 2 = 0. The second equation has the typical 
form of a transport equation and is the generalization of the classical Vlasov equation for 
charged particles with abelian interactions p^(d p — eF flv (x)dp)f(x,p) = 0, where f(x,p) 
is the classical distribution function. These two equations together give a complete de- 
scription for the Wigner operator; they are equivalent to the original field equations of 
motion. 
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2.2 Semiclassical expansion 

In order to better understand the structure of the constraint and transport equations, we 
consider their semiclassical expansion in h and their classical limit. The calculation of 
quantum corrections to this limit can then be performed in a systematic way. 

As stated in the Introduction, we will in this paper consider the electromagnetic field 
as a classical (mean) field. In this approximation the equation of motion for the Wigner 
operator and its ensemble expectation value, the Wigner function, become formally iden- 
tical. Therefore we will now work with the equations of motion for the Wigner function 
by leaving off the hats over the corresponding operators. 

The field strength F^ v in Eqs. (|16D , which has to be evaluated at the shifted argument 
x — ihsd p , is defined in terms of its Taylor expansion around x and can be expressed in 
terms of the so-called "triangle operator" A = d p -d x as 

F^x - ishd p ) = e^F^x) . (19) 

The s-integration can then be done, and expanding the result in powers of h we obtain 



sin(ftA/2) 
hA/2 



d, - eF^x) d;) + ^ 2 A 2 i^(x) d; + ■ ■ ■ , 
ehfcos(hA/2) sin(frA/2)\ „ , . ^ 

IV*, P) = P , + T (^_AJ _ _L_M j FiJx) a; 

= p„-^I\F w (x) &; + ■■■ , (20) 

where the dots indicate corrections from higher orders of h or, equivalently, of the deriva- 
tive operator A. 

Expanding the Wigner function similarly in powers of h, 

W A {x,p) = Wi°\x,p) + hW^\x,p) + h 2 wi 2 \x, P ) + • • ■ , (21) 

and inserting these expressions into the constraint and transport equations (|TT| ) and (p!8|), 
we obtain their semiclassical expansion. In the zeroth order, there is no information from 
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the transport equation, and the constraint equation reduces to the classical mass-shell 
condition, 

(p 2 -m 2 )wi°\x,p) = 0. (22) 

This equation has two elementary solutions corresponding to positive and negative ener- 
gies, and we can write 

Wl°\x,p) = W+ {0 \x, p) 6( Po - E p ) + W^ {0 \x, p) 5(p + E p ) , (23) 

where E p = +Vp 2 + m ' 2 is the classical on-shell energy. 

To next order in h, the transport equation ([18]) begins to contribute. It yields the 
Vlasov equation 

p"(0„ - eF pu {x)d;)wi°\x,p) = , (24) 
which must be evaluated with the ansatz (p3|). The constraint equation yields at order h 

(p 2 -m 2 ) W 4 (1) (x,p) = 0, (25) 
with the solution 

Wi 1] {x,p) = W+ (1) (x, p) 5(p - E p ) + W^ {1 \x, p) 5(p + E p ) . (26) 
At order h 2 , Eq. (|i~7D becomes 



{p 2 -m 2 )wi 2 \x,p) 



\(d x - eF(x)-d p ) 2 + ^AF^d;+j(x)-d; 



wf\x,p), (27) 



where j u {x) = ej d A pp u W(x,p) is the average charge current of the scalar particles in 
the ensemble which generates the mean electromagnetic field via Maxwell's equation (^). 
Eq. (0) yields at order h 2 

p»(d, - eF^(x)^)^ 4 (1) (x,p) = . (28) 

The right-hand side of (|2"7| ) represents off-shell corrections to the mass-shell condition 
due to quantum corrections. Eq. (^) is the Vlasov equation for the first order Wigner 
operator. Quantum corrections to the Vlasov equation for W\ will arise only at order 
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h 3 . The fact that quantum corrections affect the Wigner function only at second order 
in H is a specific feature of scalar particles. In spinor QED quantum effects arise in first 
order of h due to spin interactions. 



We now carry out the po-average of the Vlasov equations fl24] , |28D . This is facilitated 



by their simple form ([23|j26|) resulting from the mass-shell conditions (^,^). For the 
corresponding three-dimensional positive and negative energy Wigner functions W 3 (see 
Eq. (|])) we obtain 

d t W 3 ±(i \x, p) ± (vV)W 3 ±!i) (s, p) + e(E(x) ± v x B(x))-V p W 3 ±{i \x, p) = , (29) 

where v = p/E p , % = 0, 1, and E and B are the electric and magnetic field components 
of F^ v . We now introduce classical particle and antiparticle distribution functions, / and 
/, via 

/«(x,p) = W^\x,p), 

P\x,p) = W 3 {l \x,-p). (30) 

They satisfy the well-known three-dimensional Vlasov transport equations for scalar par- 
ticles and antiparticles moving in an external electromagnetic field: 

d t f^(x,p) + (v-V)/ (i) (x,p) + e(E(x) + v x B(x))-V p f\x, p) = 0, 

dj^(x,p) + (v-V)/«(i )P ) -e(E(x) +v x B(x)).V p f®(x,p) = 0. (31) 

2.3 Pair production in transport theory 

In this subsection, we give an exact non-perturbative solution of the scalar transport 
theory in a spatially constant external electric field. The solution describes pair production 
due to vacuum excitation. In electrodynamics, it was first studied many years ago by 



Schwinger ||16|| , who connected the probability of pair creation with the imaginary part of 
the effective action in QED. Recently it was reinvestigated as an application of the equal- 
time transport theory. For scalar electrodynamics in the Feshbach-Villars representation, 
Best and Eisenberg [jXQj] obtained from the kinetic theory the same result as Popov did 
previously in |2T| with field operator techniques. Here we consider the pair creation of 



scalar particles directly in the Klein-Gordon representation of the kinetic theory, to show 
how the energy averaging method works in a non-perturbative case and to discuss its 
difference from the Feshbach-Villars based equal-time method. 
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2.3.1 Initial value problem 

The determination of any solution of the differential equations fllTD and flT8|) needs initial 
conditions. For the pair creation problem we should search for a free vacuum solution 
as the initial condition. For E = B = 0, the energy average of the full constraint and 
transport equations result in the following three-dimensional expressions: 

(d 2 t -V 2 + AEl)w,(x,v) = 4e(x, P ), 

d t p{x,p)- V-j(^p) = 0. (32) 

Here we defined the phase-space densities of electric charge p(x, p), electric current j (a;, p), 
and energy e(x, p) by 

p(x,p) = J dp j (x,p) = e J dp p W 4 (x,p) , 
j(x,p) = J dp j(x,p) = epW 3 (x,p) , 

e( «,p) = fd M , P ) = J<i PaP lW l{x , P ). (33) 

The simplest homogeneous solution of Eq. (|32|) is 
1 



E p 



W 3 (x,p) 
p(x, p) = 
j(x,p) = ev = e 



P 



E p 



e(x,p) = E p . (34) 

Since purely magnetic fields do not produce pairs, we can restrict our attention for the 
Schwinger pair creation mechanism to electric fields. We assume for simplicity a spatially 
homogeneous but time-dependent electric field. The spatial homogeneity of the external 
fields and the initial condition ( |34"D allow to reduce the constraint and transport equations 
to ordinary differential equations in time. The three-dimensional kinetic equations are 
obtained through the energy average as 

D 2 t +AE 2 v + e -{d t E)-V p )w^p) = 4e(*,p), 

D tP (t,p) = 0, (35) 
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where 

D t = d t + eE • V p . (36) 

To obtain a closed equation of motion for the Wigner function W 3 (t, p), we eliminate 
the energy density e from the first equation of (j3l)| ). To this end we multiply the transport 
equation ([18|) by p from the left and then integrate it with respect to po : 

D t e(t,p) = ^(^E)-V p A + ^(9 t 2 E)-V p + ep-E)^ 3 (t,p). (37) 



We then apply the operator D t on the first equation ( |35| ) and combine it with Eq. (|37| ) to 
eliminate e. We obtain 

(A 3 + 4£ P 2 A + 4eE-p) W 3 (t, p) = , (38) 

which, together with the initial condition 

W 3 {t = -oo) = , (39) 
tip 

can be solved as an initial value problem for the Wigner function. We assume that 
the external electric field is switched on adiabatically in the far past and switched off 
adiabatically in the far future, E(i = — oo) = E(t = oo) = 0. 

The momentum derivative hidden in the operator D t complicates the solution of 
Eq. (|3^). We follow [|TJJ] and use the well-known method of characteristics in order to 
separate the momentum derivative and obtain an ordinary differential equation in time. 



One introduces a test Wigner function |L4| through 



W 3 (t,p) = J d 3 Po ^,Po)5 (3) (p-p(t,Po)), (40) 

where the function p(t, p ) is a solution of the classical equation of motion for a particle 
with initial momentum p in an external electric field: 

= eE(t) . (41) 
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Its explicit form reads 

p(*,Po) =Po + e f dt'E{t'). (42) 

J — oo 

Substituting fl4"0|) into fl38|), the momentum derivative is absorbed into the classical motion, 
and the partial differential equation is converted into 

J d 3 p [($ + 4E 2 p d t + 4eE-p)uj 3 (t,p )} 6® (p - p(t, p )) = . (43) 
This is solved if u 3 (t, po) satisfies the ordinary differential equation 

(df + 4E 2 p (t, po) $ + 4eE-p(t, p ))uj 3 (t, p ) = (44) 
with the initial condition 

u 3 (t = -oo,p ) = — , (45) 

where p(t, p ) satisfies Eq. ( fi2~D and E p (t, p ) is the correspondimg on-shell energy. 

Please note the time dependence of the momentum p and the particle energy E p 



arising from the classical equation of motion: when inserting the solution of (44) into (40) 
in order to construct the Wigner function W 3 from the test function u 3 , all the classical 
time dependence of W 3 resides in these functions p(t,Po), E p (t,po), while the additional 
time dependence from quantum effects is described by the differential equation ([14]) for 
u 3 . 

For later it will be useful to similarly introduce a test charge density g(t, p ), a test 
charge current j7"(t,po), and a test energy density £(t, p ); for example, 

p(t, p) = J d 3 p g(t, po) 5 (3) (p - p{t, po)) • (46) 

They satisfy for arbitrary initial momentum p the following equations: 

d t Q{t, Po) = , 

d t £(t,p ) = E-J(t,p ), (47) 



which follow from Eqs. (poj) and (|37|), respectively. Their physical meaning is evident: The 
first equation shows that the net charge density always vanishes due to the homogeneous 
initial conditions (|34|) . The second one is just Poynting's theorem of energy-momentum 
conservation. 
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2.3.2 Pair density 



As pointed out in |21j , the problem of pair creation in a homogeneous, but time-dependent 
electric field can be mapped onto the quantum mechanics problem of an oscillator with 
variable frequency. This will be exploited in the following treatment. We introduce an 
auxiliary function ((t) via 

uJ 3 (t) = \((t)\\ (48) 



in terms of which the pair production problem (^4| , ^5|) ) is reduced to the solution of a 
quantum oscillator problem: 



dfC + E 2 p (t) C = , |CC* = -oo)| = l/y/E p . (49) 

Since we are interested in the total pair production yield at time t — > oo, it is suffi- 
cient to study the asymptotic solutions of this equation. Due to its similarity with the 



time-dependent barrier-potential problem [22] in non-relativistic quantum mechanics, the 



asymptotic solutions can be easily written down in WKB approximation: 



C(t--oo) = e"^ 7VV> 

C(t^+oo) = C ie - iE ^ + C 2 e iE ^ . (50) 
Here we have already taken into account the initial condition for (, and E~°° and E°f are 



the asymptotic particle energies E 00 = Jm 2 + p§ and E°f = Jm 2 + (p(£ = oo, Po)) 2 - 



The condition for applicability of the WKB approximation [23 



d t E p = ^ « El , (51) 

is satisfied in the limit t — > ±oo since the electric field vanishes in this limit. 

From these asymptotic solutions we see that the test Wigner function in the limit 
t — > +oo contains both oscillating and non-oscillating parts: 

W3 (t -> oo) = [CC* ^ oo)| 2 = (|d| 2 + |C 2 | 2 ) + (dCZe-**?* + c.c.) . (52) 
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The created pairs are separated and accelerated by the electric field in opposite directions, 
thereby generating a current. Writing this current as J{t) = ev(t)n(t), where n(t) is the 
total particle density (positive plus negative particles), and comparing with its expression 
in terms of the test Wigner function, namely, J{t) = ep(t) ou 3 (t), we find 

n(t) = E p (t) co 3 (t) . (53) 

Taking the time average of Eq. ( p2[ ) which removes the rapidly oscillating parts, we thus 
find for the asymptotic value of the total particle density 



n(t = ^)=E^[\C 1 \ 2 + \C 2 \ 2 ). (54) 
Due to the conservation law 



(d t ()e-(dtO( =0, (55) 



which is easily derived from (fJ5j), the expression in the square bracket is the same at 
t = +oo and t = — oo and can thus be evaluated with the initial conditions ([5(]). We 
obtain 

I^i| 2 -I^| 2 = ^- (56) 
p 

Thus the oscillator is fully characterized by the ratio 



i a 



12 



12 



(57) 



ICil 

which can be interpreted (and calculated) as the transmission coefficient of the non- 
relativistic barrier-potential problem [p^| . In terms of this ratio, the asymptotic particle 
density is given as 

n(t = oo) = = 1 + 2^— . (58) 

1 — r 1 — r 

The first term on the right-hand side is a vacuum contribution and stems from the Dirac 
sea of charged particles; it can not be measured and must be removed by renormalization 



19| . The second term arises from the pair creation. The factor of 2 reminds us that both 
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particles and antiparticles contribute to the pair current. The asymptotic pair density in 
phase-space is thus finally obtained as 

r 

Upairit = 00) = y— ^ , (59) 



which is in full agreement with the results derived both from field theory |21| and from 



the transport theory in Feshbach-Villars representation [19]. 

Even though the final results of the two approaches to scalar transport theory, the 
equal-time and the energy averaging methods, are the same, we would like to point out 
one essential difference between these two procedures: The equal-time formulation has to 
rely on the Feshbach-Villars representation, since it requires field equations which contain 
only first-order time derivatives. As a result, the Wigner operator in the equal-time 
approach is a 2 x 2 matrix and first must be decomposed into its "spinor" components, 
in an analogous way as required for spinor QED which will be discussed in the following 
section. This matrix structure leads to a set of coupled kinetic equations for the "spinor" 
components of W3. The energy averaging method works directly with the covariant Klein- 
Gordon equation. This is a scalar equation, and no such complication arises. 

3 Spinor electrodynamics 

3.1 Covariant version of the BGR equation 

We now turn to the investigation of the transport theory for spin-- particles interacting 
with an electromagnetic field. We start from the lagrangian density 

C = + ieAfj) - m)^ - \F^F^ V , (60) 

which gives rise to the Dirac equations for the complex fields tp and ip, 

i^idn + ieAfj,)^ = mi/), 

i(dfj, — ieA^i/j^ = —mip . (61) 
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The spinor Wigner operator is the four-dimensional Fourier transform of the gauge in- 
variant density matrix 

I 

$4(x, y) = ij) (x + \yj exp ie / ^ ds A(x + sy)-y ijj (x — \y\ . (62) 

Unlike scalar electrodynamics, $4 is now a 4 x 4 matrix in spin space, and the Wigner 
operator W is no longer self-adjoint. It behaves under hermitian conjugation like an 
ordinary 7-matrix, 

Wl(x,p) = ~? W 4 (x, p) 7o . (63) 

The evolution equation of the covariant Wigner operator again follows from the field 
equations of motion. Calculating the first-order derivatives of the density matrix and using 
the Dirac equations in a similar way as in the scalar case, we find a pair of equations for 
$4 in terms of the operators D and II, 

(iD^x, y) - itt^x, y)) ^(s, y) = -oti7°$ 4 (x, y) , 

(±D lx (x,y)+iIl lx (x,y))$ 4 (x,y)^ = im^(x,y)^° . (64) 
The 70-matrices have been included in order to facilitate comparison with the three- 



dimensional equal-time approach of BGR |L4 



Multiplying (^) by another 7 -matrix from the right and taking the Wigner transform, 
we derive the following kinetic equations for the Wigner operator: 

{^UD^x^p) - zII M (x,p)) 7°7 M W / 4 (x,p)7° = -im^ W 4 (x,p)^ , (65) 
(jjhD^x, p) + in M (x, p)) W 4 (x, p)Y = imW 4 (x, p) . (66) 

With an eye on the semiclassical expansion below, we have again displayed the ^-depen- 
dence explicitly. We note that these two equations of motion in phase space, like the two 



Dirac equations (|61| ) in coordinate space, are adjoints of each other. Therefore, either one 



of the Eqs. fl65|,|66D provides a complete description of the Wigner operator. By adding and 
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subtracting these two equations, respectively, we can get two different, now self-adjoint 
equations: 

i^ M { 7 V, WAl^-in^Y, Wa1°] = -imfr , W 4l °], (67) 

|^[ 7 V, W47 ] - rtV[ 7 V, W47 } = -im{ 7 , ^47°} • (68) 

Obviously ( |67| ) and (^) are symmetric under an exchange of commutators and anticom- 
mutators. Their self-adjoint properties are easily proven with the aid of Eq. fl6"3|). Since 
these two self-adjoint equations are equivalent to either of the equations ( |65| , |66D , we have 
now three choices for a full description of the Wigner operator: Eq. fl6"5|), Eq. (j66|) , or 
Eq.(^) together with Eq. (|68|) . Eqs. (|67|) and (|6~8"D do not separately provide a complete 
description. 

We would like to note another important feature of Eq. (|67D : The commutator in the 
second term of the left-hand side leads to the automatic disappearance of po contained in 
the operator Uq(x,p). Since the only remaining po-dependence is in W4, this renders the 
energy average very simple. The result of an integration over po is just the BGR equation 
from the equal-time formulation, 

MW 3 (x,p) = -§M>{7°7, Wg(x,p)}-iII.[7 7, W 3 {x,p)]-im[y , W 3 (x,p)}, (69) 

with 

1 

D t (x,p) = d t + e j i dsB(x + ishV p ,t)-V p , 

"2 

D(ac,p) = V + e A' dsB(x + ishV p ,t) x V p , 

II(x, p) = p — ieh 2 ds sB(x + ishV p ,t) x V p , (70) 

"2 

where we have employed the BGR definition of the equal-time Wigner operator [14] : 

W 3 (x,p) = J dp W,(x,p) l0 (71) 
= J d 3 y e~ ip - y V (x + §y, t) exp -ie f \ ds A(x + sy, t) • p ^ (x - §y, t) . 
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Clearly equation (|67[) is the covariant version of the BGR equation. However, as we 
have stressed, (|5"7|) is not complete. Only together with the self-adjoint equation flB"8|) 
one obtains a full description for the Wigner operator. Therefore the covariant version of 
the BGR equation is only one part of the full covariant transport theory. The question 
thus arises whether on the three-dimensional level the BGR equation can be complete or 
not. To answer this question we will reduce Eq. (|68| ) to the three-dimensional level by 
performing an energy average, too, and study its implications. 

3.2 Energy average of full covariant theory 



The two self-adjoint equations (^3) and fl68|) have a complicated structure due to the 
occurrence of commutators and anticommutators. But since they are equivalent to either 
equation (|65|) or (|66"D, we can consider equation fl6"5|) instead. To simplify the calculation 



further, we even remove the 70-matrices and thus obtain the VGE [13] equation: 



'y (n M + ±Mj - m] W 4 (x, p) = . (72) 

As in the scalar case, we will from now on consider the electromagnetic field as classical 
and pass from the Wigner operator to the Wigner function by leaving off the hats. 

The Wigner function in spinor electrodynamics is a complex 4x4 matrix. VGE dis- 
cussed its spinor decomposition and derived a set of coupled equations for the components 
of W4. Because of their characteristic transformation properties under Lorentz transfor- 
mations, it is convenient to choose the 16 matrices Tj = {1, i^ 5 , 7^, 7^75, \o^v\ as the 
basis for an expansion of the Wigner function in spin space: 

W 4 (x,p) = \ [F(x,p) + t l5 P(x,p) + ^V»(x,p) + 7/i 7 5 ^(x,p) + \a^{x,p)] . (73) 

All the components F, P, V^, A^, and S^ v are real functions since the basis elements Tj 
transform under hermitian conjugation like W4 itself, r] = 701^70 (see Eq. (|63D ). They 
can thus be interpreted as physical phase-space densities. The expansion ( [73] ) decomposes 
the VGE equation into 5 coupled Lorentz covariant equations for the spinor components. 
Since these components are real and the operators D and II are self-adjoint, one can 
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separate the real and imaginary parts of these 5 complex equations to obtain 10 real 
equations: 

IF% = mF , 
TiD^A^ = 2mP , 



and 



-TiD^P + e^pWS*? = 2mA fl , 
\h{D^V u - DyVy) + e^pWA" = mS^ , 

hD»V„ = 0, 



n% = o, 

\hD,F = -U v S UfM , 



(74) 



4 1 br: ^tuijp- Ly u 



(75) 



Peforming an energy average of the expansion of W4, (|73|) , we obtain a full mapping 



from the four- to the three-dimensional components introduced in [14] 



W 3 (x,p) 



namely, 



fo(x,p) 
h(x, p) 

f2(x,p) 

fs(x, p) 
go (^, p) 
gifop) 



f (x, p) + 75/1(2:, p) - 27075/2(2, p) + 7o/3(x, p) 
+757"go(^, P) + 7o7-gi(z, p) - iT'fefo p) - TsTgsfo p) 

dp V (x,p) , 

dp A (x,p) , 
dp P(x,p) , 
dp F(x,p) , 

dp A(x,p) , 
dp V(x,p) , 
-Jd Po S 0l (x,p), 2 = 1,2,3, 
±e^ fc / dp S jk (x,p), 2 = 1,2,3. 



(76) 



(77) 
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As shown in fl4] , P0| , the physically interesting densities such as charge current, energy 
momentum tensor and angular momentum tensor, can also be expressed in terms of these 
spinor components. For example, fo and go are the charge density and the spin density, 
respectively. 

The VGE equations d7^ , [75| ) can be divided into two groups. Group I contains the 
operator U (x,p) which involves poi group II contains no po-dependence except for the 
one in the Wigner function itself. For group I the energy average is straightforward. The 



result is in complete agreement with the spinor decomposition of the BGR equation \TA 

ft(A/o + D. gl )=0, 
%{D t h + D-g ) = -2mf 2 , 
hD t f 2 + 2n-g 3 = 2mA , 
fiD t f 3 -2U-g 2 ) = 0, 

h(D t g + da) - 2n x gl = o , 

h(D tgl + D/ ) - 2n x go = -2mg 2 , 
h(D tg2 + D x g3 ) + 2II/3 = 2m gl , 

ft(Ag 3 -Dxg 2 )-2n/ 2 = 0. (78) 

Due to the additional po-dependence from the operator Ho, the equations in group 
II can not be completely reduced to expressions from the set of three-dimensional com- 
ponents given in Eq. (|77|) . They contain additionally higher p - m oments of the four- 
dimensional components, 

J dpop V - II • gl + flo/o = mf 3 , 
J dpopo A + n ■ go - fto/i = , 

J d PoPo P + lhD-g 3 + n f 2 = o, 

J dpop F - \KD ■ g 2 + n / 3 = mfo , 

J dpopo A + \KD x g! + lift - flogo = -mg 3 , 

J dp p V - \KD x go - II/o + n gi = , 

J dpoPo s % - ±;ld/ 3 + n x g3 - n g 2 = , 

J dp p Sj^ei -hBf 2 + 2Uxg 2 + 2n g 3 = 2mg , (79) 
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where we introduced the three-dimensional operator 



n (x, p) = ieh / dss E(x + isTi\/ p , t)-V p . 



(80) 



This second set of coupled equations for the components of W 3 does not include the 
operator D t . It corresponds to the energy average of the spinor decomposition of the 
second self-adjoint equation (|68"D, where D t drops out due to the commutator in the first 
term. These equations are therefore constraint equations rather than equations of motion. 
They arise in addition to the BGR equations. 

3.3 Classical limit 

Before discussing the non-perturbative consequences of Eqs. (|79|) we show how they con- 
strain the transport equations in the classical limit. As % — > 0, the original VGE equation 
can be written in the following quadratic form: 



This shows that the classical Wigner operator and hence all its spinor components are all 
on the mass shell, and thus they have positive and negative energy solutions similar to 





Eq.(H), 



ft w (x, p) 6( Po - E p ), f^>( x , p) 6( Po + E p ) , 
g+ (0) (x, p) S( Po - E p ), gr (0) (x, p) S( Po + E p ) , 



i = 



0,1,2,3. 




With this the remaining p - m t e g ra l s i n Eqs. (|79|) can be done, and everything can be 
fully expressed in terms of the three-dimensional functions (|77|). Eqs. (|79|) constribute 
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altogether six independent constraints for the positive and negative energy parts of the 
spinor components: 

,±(o) _ ^ P'gp" ' 
Ji - =t p > 

/ 2 ±(0) = o, 

f ±(o) _ , m ±( ) 

„±(°) 



±|-/„ ±(0) 



±(o) P x gj (0) 

g 2 = , 

m 

P 2 ±(0) / ±(0K 

±(o)__ u jfc P§o - IP-go JP 



gf (U> = v " ° u ^ . (83) 



mEp 



In addition, the classical limit |2(J of the BGR equations in (|78|) gives another four con- 
straints. Three of them are, however, already included in namely, the equations for 
ft^\ gi~ an d g^ ''- The last one, = p-g^°' ) /' m , is not independent of the above 
equations either, but can be obtained through a linear combination of Eqs. (|83| ) . There- 
fore, Eqs. (|83|) provide a complete set of constraint equations for the three-dimensional 
spinor components in the classical limit. There remain only two independent components, 
for example, the charge density and the spin density g^ . The classical limit of the 
BGR equations yields only the positive energy parts of the latter 4 constraints, and misses 



the first and the third additional equations from (|83|). 

The classical transport equations for /q and arise from the first order in h of the 
kinetic equations. From the BGR equations in (|78|), we have 

(d t + eE-V p )/ ±{0) + (V + eB x V p ) • gf {0) = , 
(d t + eE-V p )g 3 ±(0) - (V + eB x V p ) x g 2 ±(0) 

+ P h dt + e E-V p )A ±(0) + (V + eB x V p )-g ±{0) ) = . (84) 

Using the constraint equations, some straightforward manipulations lead to the following 
decoupled Vlasov-type equation for the charge density, 

cVo ±(0) ± v-V/ ±(0) + e(E ± v x B)-V p / ±(0) = , (85) 
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and the three-dimensional kinetic equation for the spin density, 



^g ±(0) ± (v. V)g ^ + e 



±(0) 



'E± V X B)-V r 



±(0) 



E 2 

p 



> go ±m )E - (E-p)a 



± f B x gf <»» 







*6) 



Through the energy averaging method, it is easy to prove that the above equation is 
just the three-dimensional formulation of the covariant transport equation for the classical 



axial vector A^\ 



P"(d„ - eF^d;)A^ = eF ap A<*® 



which can be derived from the first order in % of the linear equations (74) and (|75D . 
From the discussions in , this transport equation for A^' is equivalent to the covariant 
Bargmann-Michel-Telegdi (BMT) equation |5], IB, [23]] for a spinning particle in a constant 
external field, 

m— = eF^(r)s v (r), (88) 
dr 

where s M = Ajp / (A^ ' -A^) 1 / 2 is the covariant spin phase-space density. Therefore, the 
kinetic equation (j36|) for the spin density go = ^°' ) can be recognized as the three-dimensional 
BMT equation. It describes the procession of the spin-polarization in a homogeneous 
external electromagnetic field. 

Equations (^) and (|5B| ) can be further simplified by introducing classical particle and 
antiparticle distribution functions in an analogous way to the case of scalar QED: 



fi(x,p) 
gifop) 

/iO,p) 



+(0)/ \ 

f- (0) (x,-p), 

-(0)/ \ 



i9) 



The resulting equations for the particle distributions read 



dtfo + v-V/o + e(E + v x B)-V p /o = 



(90) 
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d t go+ (v-V)g + e[(E + v x B)-VJ go- ^ [(p-go)E - (E-p)g l + x g = ;(91) 

J E p l J E p 

those for the antiparticle distributions fo and go differ only by a minus sign in front of 

the electric charge e. 

The other particle and antiparticle distribution functions can be derived from the 
constraints (|53|): 



, PJto j Pgo 

ftp ftp 

/ 2 = / 2 = o, 

m _ m - 

J3 — TTJo > /3 — — "FT JO , 

P t - P f 

gi — "Frio j gi — "Frio > 

ftp -ftp 

P x go - P x 



g2 



^2„_ _ /V,.„. Z?2 r 



£pgo-(P"go)p _ ^pgo-(p-go)p fno , 

gs = " v • (92) 



E p m " E p m 



The discussion above demonstrates the incompleteness of the equal-time BGR formu- 
lation, (p9|) and (|Tg|), in the classical limit. Without the additional constraints (|79|) one 
has |2(J 4 independent components, /d°\ g ^ and g!j , and one obtains two groups of 



coupled equations, where the first one couples / to / 3 and the second one couples g 
to g^. The additional constraints arising from the complete covariant formulation allow 
to decouple from and g^ from gg . Therefore, we need only one scalar and one 
vector component to completely describe the classical behavior of the Wigner operator. 
In fact, already VGE ]13| pointed out in their covariant formulation that to any finite 
order in h the pseudoscalar P, vector and antisymmetric tensor S^ u components can 
be expressed in terms of the scalar F and axial vector components. Furthermore, the 
general relationship Tl^A^ = in ( [74]) shows that A is not an independent component 
either. Thus only 4 of the 16 components of the Wigner function are dynamically in- 
dependent. Through our mapping between the four- and three-dimensional components, 
which results from taking the energy average, we are thus able to construct transport 
equations with only 2 independent densities, the scalar mass density / 3 and the vector 
spin density g , to any order of the semiclassical expansion. 
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3.4 General kinetic equations in three-dimensional form 



In this subsection, we go beyond the semiclassical expansion and study the full quantum 
kinetic equations in three-dimensional form. This is important for the treatment of non- 
perturbative processes like pair creation in strong electric fields. Beyond the classical limit, 
the classical mass-shell condition is generally violated. The quantum Wigner function is 
no longer a 5-function in po located at the mass-shell energy, and for the elimination of 
the higher po- m oments of the covariant spinor components from Eqs. (|79|) we must follow 
a different path. Following the procedure from Section 2.3.1, we multiply the equations in 
group I of ( |7^J75| ) by p from the left and then take the energy average of these equations. 
We find 

A J dpoPo V + D J dp oPo V + If + J-gi = , 

A / dpopo A + ~D- J dp po A- I fx- Jg = 2m J dp p P , 

A / dpopo P + If 2 + 2 K-g 3 = -n J dp p e ijk S jk ei - 2m J dp p A , 

\ [A I dp p F + If 3 — 2 K-g 2 ] = -n- J dp p S m ei , 

A [dpopo A + D J dp p A - 7g - J/j + 2 Kx gl = 2 II x J dp p V , 

| [A j ' dpoPoV + D J dp po V + Jgi + J/ - 2 Kxg = 

-Ilx J dp p A + m J dp p S^e; , 
A J dpoPo S^e, - ^Dx J dp p e ijk S jk ei - Jg 2 - Jxg 3 - 2 K/ 3 = 

2U. J dpop F -2m J dp p V , 
| [| A / rfpo Po e l]k S jk e t + BxJd PoPo S^e, + Ig 3 - J x g 2 - 2 K/ 2 = 

njdpoPoP, (93) 

with 

i 
_• 



J(x, p) = ie / i c?s s (<9 4 E)(x + isV p , £)-V p , 

l l 

3(x, p) — ie ds s (<9 t B) (x + isV p , t) x V p — e ds E(x + isV p , t) , 

"2 •'""a 

K(x,p) = e [\dss 2 (d t B)(x + isV p ,t) x V p -Me ds s~E(x + isV p ,t) . (94) 

"2 "2 
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We now combine Eqs. (|79|) and (^) and eliminate all higher p - momen ts. In order to 
arrive at a set of independent contraints we also use the BGR equations and remove all 
information already contained in them. After a straightforward but tedious calculation 
we finally derive the following constraints: 

Lf + M ■ gl = , 
Lfi + M • g = , 
L/ 2 + 2N-g 3 = 0, 
L/ 3 -2N-g 2 = 0, 
Lg -M/ 1 -2Nxg 1 = 0, 
L gl + M/o + 2N x g = , 
Lg 2 + M x g 3 - 2N/ 3 = , 

Lg 3 - M x g 2 + 2N/ 2 = , (95) 

with 

L(x, p) = ie / ^ dss (v x B(x + isV p , tfj ■ V p , 

"2 

M(x, p) = ie / ' rfs sV(E(x + isV p , t)-V p ) + e / ' ds (e(x + isV p , t) - E(x)J , 

~ 2 _ 2 

1 

N(ar, p) = |e rfs (V p -V) E(x + isV p , t) + K(x, p) . (96) 

"2 

The complete set of kinetic equations for the Wigner function Ws(x,p) (resp. its spinor 
components) is given by the BGR equations ( fT^ j and the constraints (\p3[). These equations 
together form the complete transport theory in its three-dimensional formulation. They 
are fully equivalent to the covariant kinetic approach defined by either of the two equations 
( [55] ) or by the pair of self-adjoint equations (^,^). While the equations in fl9"5| ) do not 
involve time derivatives and thus do not by themselves describe transport, they provide 
essential constraints for the Wigner operator. 

The constraints (|9~5|) hold for arbitrary external fields. For homogeneous fields, the 
spatial derivatives of the electric and magnetic fields E and B in the operators L, M and 
N vanish, and the constraints (|95| ) disappear identically. This shows that the equal-time 
BGR theory is in fact complete for the case of a spatially constant external field, i.e. for 
the case studied in ||14| . The application of the BGR theory in ]14| to the problem of pair 
creation in spatially homogeneous but time-dependent electric fields is therefore safe. 
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4 Conclusions 



Relativistic kinetic theory, formulated in terms of gauge covariant Wigner operators, pro- 
vides a useful method to study transport problems in quantum theory, especially the 
phase-space evolution of high temperature and density plasmas. There exist two paths 
to formulate such a transport theory, starting from 4-dimensional or from 3-dimensional 
momentum space, respectively. In this paper we have shown a way to connect these two 
approaches directly and arrived at the 3-dimensional formulation by taking an energy 
average of the covariant 4-dimensional formulation. 

We demonstrated the method explicitly for scalar and spinor QED with external elec- 
tromagnetic fields. In scalar QED, our approach started directly from the Klein-Gordon 
equation and thus avoided the complications connected with the 2x2 matrix structure of 
the Feshbach-Villars representation ||15|| . Using the energy averaging method and taking 



the semi classical limit, we directly arrived at the classical Vlasov equation for on-shell 
particles. We pointed out that there are no first-order quantum corrections due to the 
absence of spin. This may partially explain the success of the classical Vlasov approach 
as a good approximation for the description of many quantum systems involving scalar 



fields J25]. 



Beyond the semiclassical expansion, we could still perform the energy average ana- 
lytically for the case of spatially constant electromagnetic fields, and we recovered the 
correct result for the pair creation rate in a homogeneous electric field. Our calculation 
was much easier than the one based on 3-dimensional transport theory in the Feshbach- 
Villars representation, due to the Lorentz scalar nature of the Wigner function employed 
by us. 

In spinor electrodynamics, we concentrated on the question of completeness of the 
3-dimensional (BGR) transport theory. Employing the energy averaging method, we 
identified a covariant version of the equal-time BGR equation, but found that this co- 
variant equation is not equivalent to the original Dirac equation. Integrating instead 
the covariant VGE equation, which forms a complete description for the spinor Wigner 
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function, over the energy p , we derived a set of 3- dimensional (equal-time) transport 
and constraint equations which is complete and contains the BGR equations as a subset. 
The additional equations can be understood as resulting from the energy average of the 
4-dimensional constraint or generalized mass-shell condition. They cannot be obtained 
by the direct equal-time approach of BGR. In the semiclassical limit they allow to reduce 
the number of independent spinor components of the Wigner function from 8 in the BGR 
approach to 4, one scalar charge density and 3 vector components of the spin density. 
The equations for the charge density and for the spin density decouple, and each of them 
satisfies a Vlasov-type transport equation in (6+l)-dimensional phase space. We also 
derived the additional constraints explicitly for the full quantum case and showed that 
they vanish only for homogeneous external electromagnetic fields. 

Either of the two formulations has its advantages and disadvantages. The approach 
based on 4-dimensional momentum space (8-dimensional phase space) is manifestly Lorentz 
covariant, but setting it up as an initial value problem poses certain difficulties |l4j| : since 
the covariant Wigner function is a 4-dimensional Wigner transform of the density ma- 
trix, its calculation at t — — oo requires knowledge of the fields at all times. Since this 
knowledge does not a priori exist, the covariant transport equations can only be solved 
with phenomeno logically motivated forms for the Wigner function in the far past. The 
formulation in 3-dimensional momentum space (i.e. (6+l)-dimensional phase space) does 
not have this problem: it requires the density matrix only at equal times, and the Wigner 
function at t — — oo can be directly calculated from the fields at t = — oo. Setting up the 
initial value problem is therefore straightforward in this approach. For the calculation 
of the pair creation rate in homogeneous electric fields this seems to be crucial: also our 
calculation was based on the 3-dimensional version after taking the energy average. 

In both formulations the transport equations are supplemented by constraint equa- 
tions. In the covariant approach this is essentially one 4x4 matrix equation, and it is 
easily interpreted as a generalized mass-shell condition. In the equal-time approach the 
constraints are less transparent, and their derivation via the energy averaging procedure 
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is actually rather tedious. For homogeneous external fields, however, they disappear to all 
orders in h, so for this particular case the equal-time approach has a decisive advantage 
over the covariant approach where this does not happen. In all other cases a large number 
of constraints, Eqs. (p5|) , have to be solved together with the BGR kinetic equations (ff8|). 
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